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In this paper, by using two dimensional (2D) Hubbard models with 7r-flux phase and that on a 
hexagonal lattice as examples, we explore spin-charge-separated solitons in nodal antiferromagnetic 
(AF) insulator - an AF order with massive Dirac fermionic excitations (see detail in the paper). 
We calculate fermion zero modes and induced quantum numbers on solitons (half skyrmions) in 
the continuum limit, which are similar to that in the quasi one-dimensional conductor polyacetylene 
(CH)i and that in topological band insulator. In particular, we find some novel phenomena : thanks 
to an induced staggered spin moment, a mobile half skyrmion becomes a fermionic particle; when 
a hole or an electron is added, the half skyrmion turns into a bosonic particle with charge degree 
of freedom only. Our results imply that nontrivial induced quantum number on solitons may be a 
universal feature of spin-charge separation in different systems. 
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The Fermi liquid based view of the electronic proper- 
ties has been very successful as a basis for understanding 
the physics of conventional solids. The quasi-particles of 
Fermi liquid carry both spin and charge quantum num- 
bers. However, in some cases, spin-charge separation oc- 
curs, providing a new framework for thinking about the 
given systems. It indicates that the systems have two 
independent elementary excitations, neutral spinon and 
spinless holon, respectively, as opposed to single quasi- 
particle excitation in conventional solids. 

The first example is electronic systems in one spa- 
tial dimension [![. The idea of solitons with induced 
quantum numbers starts with a beautiful result obtained 
in the context of relativistic quantum field theories by 
Jackiw and Rebbi0]. Based on this idea, spin-charge 
separated solitons had a lasting impact on condensed 
matter physics. In the long molecule-chain of trans- 
Polyacetylene, spin-charge separation can occur in term 
of soliton states Q. Due to induced fermion quantum 
numbers, the soliton may be neutral particles with spin 
1/2, or spinless with charge ±e. In two dimensional elec- 
tronic systems, spin-charge separation has been a basic 
concept in understanding doped Mott-Hubbard insula- 
tor related to high-T c cuprates [1, [f| . It is supposed that 
the particles can be liberated at low energies, with spin- 
charge separation being an upshot in the "resonating va- 
lence bond" (RVB) spin liquid state which is proposed 
by Anderson as a new state of matter Recently topo- 
logical band insulator (TBI) has attracted considerable 
attention because of their relevance to the quantum spin 
Hall effect [1,0. It is pointed out that there may exist 
spin-charge separated solitons in the presence of n flux 
with induced quantum numbers 

In this paper we focus on a special class of antiferro- 
magnetic (AF) ordered state - nodal AF insulator, and 
we will show how spin-charge separation occurs. Nodal 
AF insulator is an AF order (long range or short range) 
with massive Dirac fermionic excitations. When there is 
no AF order, fermionic excitations become nodal quasi- 
particles. There are two examples of nodal AF insulator 



in condensed matter physics - one is an AF order on 
a honeycomb lattice, the other is 7r-flux phase together 
with a nonzero Neel order parameter. Based on these ex- 
amples, our results confirm that induced quantum num- 
ber on solitons is an important feature of the spin-charge 
separation in different systems. 

Formulation— To develop a systematical formulation, 
we start with the extended Hubbard models, 
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Here cj a and electronic creation and annihila- 

tion operators. U is the on-site Coulomb repulsion, a 
are the spin-indices for electrons, /i is the chemical po- 
tential, denotes two sites on a nearest-neighbor 
link, fiji and hji are the number operators of electrons 
with up-spin and down spin. On a honeycomb lattice, 
the nearest neighbor hopping is a constant, tij = t; on 
a square lattice with 7r-fiux phase, it can be chosen as 
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function of the extended Hubbard models is written as 
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Ci tC! and Cj. CT arc Grassmann variables describing the elec- 
tronic fields. 

Firstly let us derive long wave-length effective La- 
grangian of the hopping term in the extended Hubbard 
models. Although 7r-flux phase does not break trans- 
lational symmetry, we may still divide the square lat- 
tice into two sublattices, A and B. After transform- 
ing the hopping term into momentum space, we obtain 
Ef = 2x-\/cos 2 k x + cos 2 k y . So there exist two nodal 
fermi-points at ki = (§,§), k2 = (§,—§) and the spec- 
trum of fermions becomes linear in the vicinity of the two 
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nodal points. On a honeycomb lattice, after dividing the 
lattice into two sub-lattices, A and B, the dispersion is 
obtained in Ref . [H, [lfl, [13] . There also exist two nodal 
points, ki = 71(1,73) and k 2 = ^=(-l,--L) and the 
spectrum of fermions becomes linear near ki 2- In the 
continuum limit, the Dirac-like effective Lagrangian de- 
scribes the low energy fermionic modes for both cases 



(3) 



t x , t v , t z are Pauli matrices. We have 



where ipi = Vi7o = ( ipUA, ^ib, ipliA, tpiiB ) and 

i'2 = IpllO = ( 4>]2B, 4>]2A, 4>i2B, $12A ) [3, [3, E3] • 

7 M is defined as 70 = a (g) r 2 , 71 = a ® r y , 72 = (T ® r x , 

/l 0' 
ao = { 1 
set the Fermi velocity to be unit vf = 1. 

In the strongly coupling limit, U » t, there always 
exists an AF order in the extended Hubbard models. In- 
troducing Stratonovich-Hubbard fields for the spin de- 
grees of freedom [l3[ , we obtain the partition function as 
- P 8 d L 

Z = J VcVcVBe Jo , where the Lag rangian is given 
by 
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with Pauli matrices a =(a x , er a , er z ). Here B_, is a vector 
denoting spin configurations, Bj = |_Bj| tij where \Bj\ = 
(j)o represents the value of localized spin moments and rij 
is a unit vector describing the Neel order parameter. In 
AF ordered state, the mass gap of the electrons is given 
as m = (f>all. 

Then starting from Eq. ^ , we get the same long wave- 
length effective model of nodal AF insulator [l8, TfJ III 



£ off = 22 i^al^d^a + m(tpin ■ aipi - -^n • atp 2 )- (5) 



a = 1, 2 labels the two Fermi points 

Zero modes on half skyrmions - In this section we will 
study the properties of topological solitons. Instead of 
considering topological solitons with integer topological 
charge (skyrmions), we focus on solitons with a half topo- 
logical charge, J -^n • d x ii x d y n cPr = ±i with r = (x, 
y). Such soliton is called a half skyrmion (meron or anti- 
meron). A meron with a narrow core size (the lattice 
size a) is characterized by n = r/ | r |, r 2 = x 2 + y 2 



[HI, [H, [H, 0, [H, [3, [23, El]- To stabilize such a soli- 
tons, one may add a small easy-plane anisotropy of the 
Neel order. 

Around a meron configuration, the fermionic operators 
are expanded as 
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where 6 a k and are operators of k 7^ modes that are 
irrelevant to the soliton states discussed below. -0^ k (r) = 

( V'TaAk' V'TaBk: V'j'aAk. ^laBk ) ai'C the functions of 

zero modes. a° are annihilation operators of zero modes. 

To obtain the zero modes, we write down two Dirac 
equations from Eq.© 



and 



idx'Yiipi + id y 72^1 + mn ■ o-tp° = 



Firstly we solve the Dirac equation for i/i®. 

6(5) 
6(5) 

V 6(5K 



(7) 
(8) 

With the 



ansatz ip® 



we have 



9x6 = 6, 3z6 = 6, 
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9x6 = — r- + 6, #56 = 



where r = | r | (cos 9, sin 9) and i 
has been obtained in Ref. [29fl as 



h6 

The solution 



6(5) = 6(5) = 0, 6(5) - 6(5) = x*Ki(x) (10) 
where Ki(x) is the modified Bessel function. So the 



solution of ipi becomes 



i 2 if 1 (x) 

Vo 



To solve V27 we transform the equation idijiip® 
mn • cr-02 = into 



mUn ■ o-^lf- 1 = 0, 



where U 



i7T7o/2 
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U^iU- 1 = -7< and U^^U = 

\ 



Then the solution of i s obtained as 
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It is noticable that from above solutions of zero modes, 
the components if^iA, an( l ^123 are zero. 

Topological mechanism of spin-charge separation - For 
the solutions of zero modes, there are four zero-energy 
soliton states | sol) around a half skyrmion which are 
denoted by | 1+)® | 2+), | l_)ig | 2_), | 1_)® | 2+) 
and I 1+)® I 2_). | 1_) and | 2_) are empty states of the 
zero modes ip® (r) and , 02( r ); I 1+) an d | 2 + ) arc occupied 
states of them. Thus we have the relationship between 
a° a and | sol) as 

a? I 1+) =| 1_), a? I 1_> = 0, a° 2 \ 2+) =| 2_), a° 2 \ 2_) = 0, 

(12) 
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Firstly we define total induced fermion number opera- 
tors of the soliton states, N F = E N a ,F with 



Na.F = 



: 0«* : d 2 r (13) 

= (a^) f a°+X;(^ak-4Lk*«k)-^ 

: ?/4V> a : means normal product of V^V^. From the re- 
lation between aP a and | sol) in Eg. (p~2|) . we find that 
| 1±) or | 2±) have eigenvalues of ±^ of the total induced 
fermion number operator N F , 



N hF \l±) 
No 



±i|l±), AT 1)F |2 ± ) = 0, 
f 2 , F \2±) = ±l\2 ± ), N 2 . F \1 ± )=0. 
id 

tor is staggered spin operator, 



(14) 



Another important induced quantum number opera- 

I 



-^{iB^Tls) + (^{ 2j4 -0 T 2A + V>| 2 B^12B - $;2A^i2A - 

■0| 2B ^|2b)] : rf 2 r. For the four degenerate zero modes, it 

can be simplified into S? v ^ | sol) = ^(N2 tF — ^Vl,f) | sol). 
Let us show the detailed calculations. From the zero so- 
lutions of V^iai VqiB> V'^A an< ^ V'^B' we obtain four 
equations 



(15) 
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Using above four equations, we obtain 
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Then we calculate two induced quantum numbers de- 
fined above. Without doping, the soliton states of a half 
skyrmion are denoted by | 1_)® | 2 + ) and | l+)<& | 2_). 
One can easily check that the total induced fermion num- 
ber on the solitons is zero from the cancelation effect 
between two nodals N F | l_)(g) | 2+) = N F | 1+}® | 
2_) =0. It is consistent to the earlier results that forbid 
a Hopf term for the low energy theory of two dimen- 
sional Heisenberg model[30|]. On the other hand, there 
exists an induced staggered spin moment on the soliton 
states | l_)(g> | 2+) and | 1+}® | 2_), 
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| 1_)® | 2+) = - | 1-}® | 2+), 
| 1+)® | 2_) = -\ | 1+)® | 2_). 



(16) 



The induced staggered spin moment may be straightfor- 
wardly obtained by combining the definition of S? n ^ and 
Eq.KHJ) together. 

When half skyrmions become mobile, their quantum 
statistics becomes important. Let us examine the statis- 
tics of a half skyrmion with an induced staggered spin 
moment. In CP(1) representation of n, a "bosonic 



spinon" is introduced by n = zaz with z = 



and zz = 1. Since each "bosonic spinon" z carries ±i 
staggered spin moment, an induced staggered spin mo- 
ment corresponds to a trapped "bosonic spinon" z. On 
the other hand, a half skyrmion can be regarded as 
a 7T— flux of the "bosonic spinon", J • d x n x d y ii 



d2r = h I e^^a^r 



±2- with a M 



r(z<9„ 



<9 M zz) 



To be more explicit, moving a "bosonic spinon" z around 
a half skyrmion generates a Berry phase <fi to z — > z' = 

T I where <j> = J e lll/ d ll a 1/ d 2 r = ±n. As a result, 

a "bosonic spinon" z and a half skyrmion (mcron or an- 
timeron ) share mutual semion statistics. Binding the 
trapped " bosonic spinon" , a mobile half skyrmion be- 
comes a fermionic particle. We may use the operator 
f a to describe such neutral fermionic particle with half 
spin. The relation between the zero energy states and 
the fermionic states is given as | 1+}® | 2_) = fl | 
0)/ and | 1_)® | 2+) = /j | 0}/ (The state |0)/ is 

defined through / T |0)/ = /JO)/ = 0). We call such 
neutral object (fermion with ±i spin degree freedom) 
a (fermionic) "spinon". 

Now we go away from half filling. It is known that 
when a hole (electron) is doped, it is equivalence to re- 
moving (adding) an electron. Without considering the 
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TABLE I: Quantum numbers of the degenerate soliton states. 



existence of half skyrmions, the hole (electron) will be 
doped into the lower (upper) Hubbard band. The ex- 
istence of zero modes on half skyrmions leads to the 
appearance of bound levels in the middle of the Mott- 
Hubbard gap 22| . The hole (electron) will be doped onto 
the bound states on the half skyrmion and then one of 
the zero modes is occupied. When one hole is doped, the 
soliton state is denoted by |1_) ® |2_). One can easily 
check the result by calculating its induce quantum num- 
bers. On the one hand, there is no induced staggered spin 
moment, ^ | 1_)<X> | 2_) = 0. On the other hand, 

the total fermion number is not zero, 7Vp|l_) <g) |2_) = 
— 11_) <E> |2_). These results mean that such soliton state 
is a spinless " holon" with positive charge degrees of free- 
dom. After binding a fcrmionic hole, the soliton state 
(holon) does not have an induced staggered spin moment. 
Thus the holon obeys bosonic statistics and becomes a 
charged bosonic particles. When one electron is doped, 

_i_). The induced 



the soliton state is denoted by |1+) 
quantum numbers of it are Np\l+) £* 



|2+) = +|l+)®|2+) 



and Sf 



0. Such soliton state is also a 



>)r) i 1+>®J 2+) 
bosonic particle with a negative charge but without spin 
degrees of freedom. We call such a soliton state an "elec- 
ton" to mark difference with the word " electron" . 

Finally we get a topological mechanism of spin- charge 
separation in nodal AF insulators. There exist two types 
of topological objects - one is the fermionic spinon, the 
other is the bosonic holon ( or the bosonic electon ). 

In ID system, real spin-charge separation may occur. 
As far as the low energy physics is concerned, the spin 
and charge dynamics are completely decoupled from each 
other. In 2D, real spin-charge separation in a nodal AF 
insulators can not occur in long range AF order. In the 
future we will study the deconfincment condition of spin- 
charge separated solitons and explore the properties of 
deconfincd phases with real spin-charge separation. 

Summary - By using 2D 7r-flux phase Hubbard model 
and the Hubbard model on a honeycomb lattice as ex- 
amples, we explore spin-charge separation in nodal AF 
insulator. The crus crux of the matter in this paper is 
the discovery of induced staggered spin moment S? n ^ 
on half skyrmions in nodal AF insulators. Based on such 
nontrivial induced quantum number, we classify four de- 
generate soliton states with zero energy - two of them 
(| 1_)35 | 2 + ) and | 1+)® | 2_)) represent the up-spin 
and down-spin states for a fermionic " spinon" , another 
state (|l_)(g>|2_)) represents a "holon" and the last one 
(|1+) ® |2+)) denotes an "electon". 
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